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Dark matter detection rates for supergravity models with R parity where supersymmetry is broken at a scale 
~ Mg are discussed. Non-universal soft breaking masses in both the Higgs and squark sectors are considered, and 
it is seen that these can effect rates by a factor of 10 - 100 when m^o ~ 65 GeV (xi = lightest neutralino) but 
otherwise make relatively small corrections. The b — > s + 7 branching ratio is seen to correlate with detector event 
rates, large (small) branching ratios corresponding to small (large) event rates. Effects of precision determinations 
of cosmological parameters or event rate predictions by future satellite experiments are discussed for the ACDM 
and the i/CDM models 



1. INTRODUCTION 

The dark matter problem is a particularly in- 
teresting one to examine in supersymmetry since 
models with R-parity invariance automatically 
have a dark matter candidate, the lightest su- 
persymmetric particle (LSP). Over most of the 
parameter space in supergravity models, the LSP 
is the lightest neutralino (x?), and so predictions 
for event rates for detection can be made. 

We consider here models based on supergrav- 
ity where SUSY is broken in a hidden sector 
at a scale ~ M G (where M G = 2 x 10 16 Gev 
is the GUT scale) with this breaking communi- 
cated to the physical sector by gravity [1], and 
with radiative breaking of SU(2) x U(l) occur- 
ring at the electroweak scale ~ Mz [2]. The 
simplest such model is the minimal one [MSGM] 
with universal soft breaking at Mq [1,3]. Such 
models depend on only four additional parame- 
ters and one sign. These can be taken to be Too 
(the universal scalar soft breaking mass), TO1/2 
(the universal scalar soft breaking mass), Ao (the 
universal cubic soft breaking parameter), B (the 
quadratic soft breaking parameter) and the sign 
of no (where /io is the Higgs mixing parameter 
in the superpotential, W M = fio HiH 2 ). After 



SU(2) x U(l) breaking, one may alternately re- 
place B by tan/3 = {H 2 ) / {H2) (#1,2 are the two 
Higgs doublets of SUSY), TO1/2 by the gluino mass 
(nig = {a 3 (Mz)/a G )m 1 / 2 , where a G = 1/24 is 
the GUT coupling constant) and A by A t (the 
t-quark A parameter at the electroweak scale). 

Until recently, almost all calculations on detec- 
tion of dark matter (DM) has been done within 
this universal framework. However, one expects 
the possibility of non-universalities arising from 
the following sources: (1) Kahler potential inter- 
actions can give rise to non-universal soft break- 
ing. (2) Even if universality holds at a higher 
scale e.g. the string scale M str , the running of 
the renormalization group equations (RGE) will 
generate non-universalities at Mq. (3) In break- 
ing higher rank grand unified groups to the Stan- 
dard Model (SM) group, the D terms can generate 
non- universalities . 

The number of parameters needed to describe 
the non-universality expected at Mq depends 
upon the gauge group. For example, for SU(5), 
neglecting phases and splittings in the first two 
generations (to supress flavor changing neutral 
currents) one needs nine additional observable pa- 
rameters, of which only four enter significantly 
for large segments of the parameter space. (This 
might be compared with the « 30 parameters 
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of the MSSM). Over the past two years effort 
has been made to explore this larger parameter 
space [4,5], and we will describe here some results 
that have been obtained. We will see that non- 
universal masses can occur both in the Higgs and 
sfermions sectors, and these two sectors can effect 
each other either constructively or destructively. 

There are several phenomena which effect pre- 
dicted dark matter rates, and we summarize these 
now: 

f. t-quark mass (m t = 175 GeV). The heavy 
top drives the lightest top squark, t\ , tachy- 
onic for negative A t , eliminating this part of 
the parameter space unless A t /m ~ -0.5. 

2. b — > ,s+7 decay. The current CLEO branch- 
ing ratio is B[B -> X sl ] = (2.32± 0.67) x 
10~ 4 [6] which can be compared with the 
SM calculation (with NLO corrections) of 
B[b sj] = (3.28 ± 0.33) x 10~ 4 [7]. ft 
is clear that already the experimental value 
constrains any new physics corrections that 
raise the theoretical value, and one finds 
that most of the A t j \i < region is already 
eliminated at the 95% C.L. Combined with 
the t-quark effect above, most of the p, < 
part of the parameter space has been elim- 
inated [5]. 

3. Amount of cold dark matter (CDM). The 
astronomical determinations of the various 
cosmological parameters at present have 
large uncertainties. Thus the Hubble con- 
stant, H = h (100 km/s Mpc, has the 
range of 0.5 - h - 0.75. The density of 
matter pi of type i can be measured by 
Cli = pi/ p c where the critial density p c is 
given by p c = 3H 2 /8irG N = 1.88 x l()- 29 h 2 
gm/cm 3 . The determinations for cold dark 
matter (non-relativistic matter at the time 
of galaxy formation) are in the range 0.3 
~ CIcdm ~ 0.75. Thus one has 

0.1 < n C D M h 2 < 0.4 (1) 

We will assume here that this CDM are the relic 
ncutralinos Xi- Eq. (1) represents the cosmo- 
logical abundance of x?- Terrestial detectors can 



observe the CDM in the Milky Way impinging 
on the solar system. This local density of Xi has 
uncertainties due to modeling of the halo of the 
Galaxy and the amount of machos in the Galaxy. 
We will assume here that the local density is p Xl 
= 0.3 GeV/cm 3 , though this number could be in 
error by a factor of 2 or more. 

The calculation of predicted event rates then 
proceeds as follows: one first calculates the relic 
density given by [8] 

Qoh 2 - 2.45 x 10- n f^Vf 

x N} /2 /J( Xf ) (2) 

where J(xf)—J^ f dx < av > GeV -2 , Xf — 
kTf/m x a, Tf is the freezeout temperature, iV/- 
the number of degrees of freedom at freezeout, a 
is the x? annihilation cross section, v is the rel- 
ative velocity, and <> means thermal average. 
We restrict the SUSY parameter space such that 
Q x oh 2 falls within the allowed window of Eq. (1), 
and also that the SUSY bounds from LEP, Teva- 
tron and CLEO be obeyed. Within this restricted 
parameter space, one then calculates the detector 
event rate R given by [9] 

R = (Rsi + R S D)(p x o/0.3GeVcm- 3 ) 

is events 

w yO /320fcms- 1 (3 

v Xl ' ' kg d 

where R SI = 16m x oM^M^[M N + m x o]~ 2 | 
A S i | 2 and R SD = l6m x oM N [M N + 
m x o]- 2 \ 2 J(J + 1) | A SD | 2 . Where M N is the 
target nuclear mass and J is its spin, and Asi, 
Asd are the spin independent, spin dependent 
scattering amplitudes. Note that for large Mjv 
one has Rsi ~ M^, while Rsd ~ 1/M/v thus fa- 
voring the spin independent scattering with heavy 
nuclei. 

2. p 2 Dependence on Soft Breaking Masses 

We will assume here that the first two gener- 
ations of sfermions are degenerate (to suppress 
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flavor changing neutral currents) and allow for 
non-universal soft breaking masses in the Higgs 
and third generation sfermions. We parameterize 
these at Mq as follows: 



m H 2 = 771q(1 + Si); m H 2 = ml(l + S 2 ) 



m d R 



+ 6 3 ); m. 
rog(l + <f B ) 
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ml(l + <5 4 ); 



ml(l + 5 e ); m\ L = m 2 1 (l + £ 7 ) 



(4) 

(5) 
(6) 



Here — arc the doublet of squarks, 

£l = (^l,cl) the doublet of sleptons and the ref- 
erence mass mo is taken to be the common mass 
of the first two generations. In addition there 
are the t, b and r cubic soft breaking parameters 
A t, j4o6, A ot . For grand unified models with 
GUT groups containing an SU(5) subgroup (e.g. 
SU(N), N > 5; SO(N), N > 10, E 6 etc.) and 
with matter in the usual 10 + 5 of SU(5), one 
has 



63 = 54 = 65; S e = S 7 ; A ob 



(7) 



In the following, we will limit our parameter 
space so that m n , m g < 1 TeV, tan/3 < 25 and 
-1 < Si < 1. For tan/3 in this domain, results are 
generally insensitive to S$, <5g, ^7 A Q b and A OT , 
and we will set these parameters to zero. The 
radiative breaking of SU(2) x U(l) determines /i 2 
to be 



11 2 - n 2 t 2 

2 Ml H-2 1 



t 2 -l 



-M\; t = tan/3 



(8) 



where /if = m 2 H . + and m 2 H . are the running 
Higgs masses, and E, are loop corrections. One 
finds [5] 



M 2 = 



t 2 -l 
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Figure 1. Maximum and minimum event rates 
for a xenon detector for universal soft breaking 
(solid), 62 = — 1 = —5\ (dotted) and 62 = 1 = 
-82 (dashed) [5]. 



2 z 22 1 2 - 1 1 a G ' 
+ loop corrections (9) 

where D Q = 1 - (m t /200sin/3) 2 , A R ^ A t - 0.61 
mg, = TrYm 2 , Y is the hypercharge and m 2 
in 5*0 are the scalar masses at Mq- D vanishes at 
the t-quark Landau pole, and in general is small 
(D ~ 0.25). A R is the residue at the Landau 
pole. One sees several features of Eq. (9). For 

t 2 » 1 (i.e. t~ 3) 81 does not enter sensitively 
in /j, 2 . Further, since Do 1S small, we see that £3 
and (*>4 acts oppositely to S2, i.e. the squark non- 
universalities can contribute constructively or dis- 
tructively to the Higgs non-universality. One thus 
cannot consider only one type of non-universality. 

In general, one can expect important modifi- 
cations due to non-universal soft breaking to oc- 
cur when, for some reason, the universal terms 
are small. This can occur if D = |, or when 
the residue at the Landau pole vanishes (A t = 
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0. 6m g ), or if m§ is small (i.e. if m x o is small since 
for much of the parameter space /i 2 /M§ >> 1 
and one has the scaling ~ (\/7)m,g [10]). 
These effects are enhanced for small tan /3. We 
illustrate some of these effects for DM detection 
event rates, /z 2 gives rise to small event rates and 
small [i 2 plays a key role here in that it governs 
the interference between the Higgsino and gaug- 
ino parts of the \i m the SI part of the Xi ~ 
nuclear scattering cross section. In general, large 
/j, 2 gives rise to small event rates and small /i 2 
to large event rates [11]. Fig. 1 shows the max- 
imum and minimum event rates for a xenon de- 
tector for universal and non-universal soft break- 
ing masses. Here 83 = = 84. One sees the 
non-universal effects are small for large neutralino 
masses (m^ ~ 60 GcV). For <5 2 = —1 = —Si, 
there can be a reduction of a factor ~ 10 - 100 in 
the minimum event rates (where tan/3 is small) for 
small m x a, since then by Eq. (9) ^ 2 is increased 
by the non- universalities. Correspondingly for 
<5 2 = 1 = —82, n 2 is decreased and event rates 
can be increased by a factor ~ 10. Fig. 2 shows 
the corresponding curves for Si = S 2 = S3 = 0. 
One sees that the Si — +1 dotted curve, resem- 
bles the 82 = — 1 curve of Fig. 1, and the 84 = — 1 
dashed curve resembles the <5 2 = +1 curve of Fig. 

1, as one would expect from Eq. (9). 

Fig. 3 exhibits the fact that when S3 = 1 = 84 (as 
expected in GUT models) and 82 = 1, the squark 
and Higgs non-universal effects act coherrently to 
significantly increase the maximum event rates up 
to ~ 10 event/kg d (which is the current level of 
dark matter detector sensitivity [4]). However, 
for 82 = — 1 = — <5i, the two effects mostly cancel, 
yielding event rates close to predictions of univer- 
sal soft breaking. 

3. Cosmological Parameters 

As discussed in Sec. 1, the various cosmologi- 
cal parameters are at present not well determined. 
However, future satellite experiments, MAP and 
PLANCK [12] will be able to measure the angular 
power spectrum quite accurately, and this will al- 
low the determination of the Hubble constant, the 
amount of dark matter, the cosmological constant 
etc. at the level of (1-10)% [13,14]. Such deter- 
minations would considerably restrict the SUSY 
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Figure 2. Maximum and minimum event rates 
for a xenon detector for universal softbreaking 
(solid), 84 = +1 (dotted), and 84 = —1 (dashed) 
[5]. 
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Figure 3. Maximum and minimum event rates 
for a xenon detector for universal soft breaking 
(solid), <5 3 = 1 = 84, 82 = —1 = —Si (dotted), 
and S3 = 1 = 84, S 2 = 1 = — 82 (dashed) [5]. 
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paramter space, and hence sharpen significantly 
the predictions of DM detection rates. To illus- 
trate what might be expected from these deter- 
minations, we consider two cosmological models. 

(i) ACDM Model 

One assumes there that CDM and baryonic 
dark matter (B) exists with a cosmological con- 
stant A such that the universe is flat: Vic dm + 
Qb + = 1- As an example, we assume that 
the measured central values of the parameters are 
ttcDM = 0.40, tt B = 0.05, Q A = 0.55 and h = 
0.62, which are consistent with current astronom- 
ical measurements. Then using the estimated ac- 
curacy that could be achieved by the PLANCK 
satellite [13], one finds 

VcDMh 2 = 0.154 ±0.017 (10) 

This window is much narrower then what is cur- 
rently assumed, i.e. Eq. (1). Eq. (10) produces 
two interesting results for the event rates: the 
minimum event rates are significantly raised for 
mjj ~ 60 GeV, and the upper bound on ^^h 2 
produces an upper bound on allowed values of 
m° 1 . One finds at the ler (2cr) ranges that mj 1 < 
70 (77) GeV, and by scaling this produces a corre- 
sponding bound rrig < 520 (560) GeV. In addition 
m ± ~ 150 GeV and for the light Higgs one has 

rrih ~ 120 GeV. It is interesting to compare these 
results with the reach of the upgraded Tevatron 
with 25 fb^ 1 of data where the gluino would be 
observable if m g ~ 450 GeV [15], the chargino if 
m-± ~ 235 GeV for about 2/3 of the parameter 

space and the Higgs if mj, ~ 120 GeV [16]. 

As discussed in Sec. 1, the b — > s + 7 decay 
branching ratio plays an important role in limit- 
ing the SUSY parameter space, as there is already 
some strain between the experimental branching 
ratio and the SM prediction. In addition, there 
exists an interesting correlation between the DM 
event rate R and B(b — ► s + 7). Large R occurs 
mainly for large tan/3 where destructive interfer- 
ence between the SM and SUSY contributions to 
the b — ► s + 7 occurs. Small R occurs mainly at 
small tan/3 where there is constructive interfer- 
ence in the b — > s + 7 decay and hence correlates 
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Figure 4. Scatter plot of R vs. B(b — > s + 7) for a 
Xc detector for the ACDM model with universal 
soft breaking masses and \i > 0. 



with large B (b — > s + 7). This can be seen in 
Fig. 4 where a scatter plot for R vs. B(6 — > s + 7) 
over the SUSY parameter space is shown for the 
ACDM model (assuming universal soft breaking 
masses). One sees for almost all parameter points 
that, if B(b -> s + 7) > 3 x 10~ 4 (the 1 std. lower 
bound of the SM), then R < 0.1 events/kg d, 
while if B(6 -> s + 7) < 3 x 10~ 4 , then R > 
0.05 events /kg d. Thus new CLEO data on the 
b — > s + 7 decay will have significant effects on 
predictions of dark matter event rates. 

(ii) i/CDM Model 

If neutralinos have mass of order of a few eV, 
they could represent a hot dark matter compo- 
nent to the dark matter. As an example of such 
a model, we assume that the measured central 
values of the cosmological parameters are 

U v = 0.20; VLcdm = 0.75, 
ft B =0.05, h = 0.Q2 (11) 
The PLANCK satellite would then determine 
QcDMh 2 with the following accuracy [13,14]: 
Q CDM h 2 = 0.288 + 0.013 (12) 
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For this case, the narrowing of the flh 2 window, 
narrows the separation between maximum and 
minimum event rates in the region m Xl ~ 65 
GeV, and for the lstd range (0.158 < tt Xl h 2 < 
0.301) produces forbidden gaps when m Xl > 65 
GeV (which, however, get filled in for the 2 std 
range). Since Q Xl h 2 is larger for this case then in 
the ACDM model, the upper bounds on the gaug- 
ino masses are larger. One finds for the 1 std (2 
std) bounds that m x o < 95 (100) GeV, m g < 700 

(720) GeV and for the chargino m x ± ~ 200 GeV. 
Thus for this model, the LHC would most likely 
be needed to discover SUSY. 

The above two examples are meant to be il- 
lustrative of what future astronomical measure- 
ments of the basic cosmological parameters will 
be able to achieve. In particular one sees how 
astronomical measurements can impact on accel- 
erator searches for SUSY particles. 

4. Conclusions 

We have considered here dark matter detection 
rates within the framework of supergravity grand 
unification with R-parity where SUSY breaks at 
a scale ~ Mg, the breaking being communicated 
to the physical sector by gravity. Such models 
automatically imply the existance of cold dark 
matter, and over a large amount of the parame- 
ter space in amounts consistent with astronomical 
measurements. 

The detection rates expected for cases of non- 
universal soft breaking masses have been com- 
pared with the minimal universal SUGRA mod- 
els. Non-universal effects can increase or de- 
crease event rates (depending on the signs of the 
non-universalities) by factors of ~ 10 - 100 in 
the domain m xi ^65 GeV (equivalently when 

rrig ~ 400 GeV) but generally have small effects 
at higher masses. One must also consider the 
possibilities of both Higgs and squark mass non- 
universalities, as both these enter in the analy- 
sis with comparable size. Thus the Higgs and 
squark non-universalities can either cancel or en- 
hance each other in the detection rates. 

Future b — > s + 7 decay data may play an im- 
portant role in uncovering new physics, and in- 
deed may be the first place that new physics is 



seen. Already the current data strongly restricts 
the SUSY parameter space. Thus combined the 
fact that the top quark is heavy (m t = 175 GeV) 
the current branching ratio for b — > s + 7 for- 
bids most of the ix < part of the parameter 
space, eliminating most of the high event rate re- 
gion for \i < 0. The event rates for \i < are 
then ^100 times smaller than for ix > 0, and this 
phenomena holds with or without non-universal 
soft breaking masses. As error flags go down, the 
b — ► s+7 decay may more strongly restrict the al- 
lowed regions of SUSY parameter space, making 
supergravity predictions more precise. 

One of the striking features of recent years is 
the "astro-particle connection" where develop- 
ments in particle theory effect astronomical and 
cosmological theory, and astronomical measure- 
ments can influence what is expected at high en- 
ergy accelerators. One may expect this interac- 
tion to strengthen in the future. Thus new satel- 
lites (MAP and PLANCK), balloon and ground 
based experiments should be able to measure the 
basic cosmological parameters very accurately (at 
the few percent level), and hence determine the 
amounts of different types of dark matter. This 
will put further restrictions on the allowed SUSY 
parameter space. In this connection we have con- 
sidered two examples of cosmological models, the 
ACDM model (which limit the gaugino masses to 
be m v o ~ 75 GeV, mj ~ 550 GeV and m ± ~ 150 

GeV, and the i/CDM model (where m x o ~ 100 

GeV, m 5 ~ 700 GeV and m ± ~ 200 GeV. In the 
first model the upgraded Tevatron might be able 
to see SUSY, while in the second model one likely 
needs the LHC (or NLC). In either case the astro- 
nomical measurements would correlate with ac- 
celerator searches for SUSY. There is also a corre- 
lation betweeen the B(6 — > s + 7) branching ratio 
and expected dark matter detection rates. Thus 
for B > 3 x 10~ 4 (i.e. B in range expected from 
the SM) the dark matter detector event rates are 
generally small i.e. R ~ (10 _1 — 10~ 4 ) events/kg 
d, while if B < 3 x 10~ 4 (smaller than the SM 
prediction) event rates in general will be large 
i.e. R ~ (0.05 — 1) events/kg d. Thus accelerator 
measurements would correlate with astronomical 
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searches for dark matter. 
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